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CS| ■ Abstract 

CN ■ 

The standard picture of cosmology assumes that a phase transition, associated with chiral 
symmetry breaking after the electroweak transition, has occurred at approximately 10/j, seconds 
after the Big Bang to convert a plasma of free quarks and gluons into hadrons. In this paper, 
_ , we consider the quark-hadron phase transition in a DGP brane world scenario within an effective 

model of QCD. We study the evolution of the physical quantities relevant to the physical description 
(~| I of the early universe, namely, the energy density, temperature and the scale factor before, during, 

and after the phase transition. Also, due to the high energy density in the early universe, we 
consider the quadratic energy density term that appears in the Friedmann equation. In DGP 
^ ' brane models such a term corresponds to the negative branch (e — —1) of the Friedmann equation 

■"sj" I when the Hubble radius is much smaller than the crossover length in AD and 5D regimes. We 

' show that for different values of cosmological constant on the barane. A, phase transition occurs 

and results in decreasing the effective temperature of the quark-gluon plasma and of the hadronic 
fluid. Finally, we discuss the case where the universe evolves through a mixed phase with an 
^'j^ I initially small supercooling period and monotonically growing hadronic bubbles. 

'. 1 Introduction 
> . 

^ , Over the past decade the possibility of the observable universe being a brane-world [1] embedded 

^1 in a higher dimensional space-time has generated a great amount of interest. This was motivated 

by the fact that there is a strongly coupled sector of Eg x heterotic string theory which can be 
described by a field theory living in 11-dimensional space-time [2]. The 11-dimensional world consists 
of two 10-dimensional hypersurfaces embedded on the fixed points of an orbifold and the matter 
fields are assumed to be confined and live on these hypersurfaces which are known to be 9-branes. 
After compactification of the 11-dimensional theory on a Calabi-Yau 3-fold, we obtain an effective 
5-dimensional theory [3] which has the structure of two 3-branes located on the orbifold boundaries. 
This scenario has motivated intense efforts to understand the case where the bulk is a 5-dimensional 
anti de-Sitter space. In this setup, gravitons are allowed to penetrate into the bulk but are localized 
on and around the brane [4] . It was then shown that in a background of a non-factorizable geometry 
an exponential warp factor emerges which multiplies the Poincare invariant 3+1 dimensions in the 
metric. The model consists of two 3-branes situated along the 5th dimension, compactified on a 



(N 



(N 



*email: k-atazadeh@sbu.ac.ir 
^masoud.ghezelbashQusask.ca 
*email: hr-sepangi@sbu.ac.ir 



1 



/Z2 orbifold symmetry where the two branes must have opposite tensions. The evolution equation 
followed from such a brane scenario differs from that of the standard four dimensional evolution 
equation when no branes are present [5] . The existence of branes and the requirement that matter 
fields should be localized on the brane lead to a non-conventional cosmology, necessitating a more 
concerted study. A large number of studies have been devoted to the effective gravity induced on 
the brane [6] and, in particular, a great amount of interest has been generated regarding inflationary 
cosmology [7] . Not surprisingly, the problem of the cosmological constant has become a focal point 
in the brane-world studies where, for example, in [8, 9, 10] a five dimensional action with a scalar 
field is non-minimally coupled to five dimensional gravity and to the four dimensional brane tension. 
There has also been some discussion on the localization of gravity [11]. A feature common to these 
type of models is that they predict deviations from the usual AD gravity at short distances. 

A somewhat diff'erent approach within the brane-world framework is the model proposed by Dvali, 
Gabadadze and Porrati (DGP) [12, 13]. It predicts deviations from the standard AD gravity over 
large distances. The transition between four and higher-dimensional gravitational potentials in DGP 
model arises because of the presence of both the brane and bulk Einstein terms in the action. The 
Friedmann-like equations governing the cosmological evolution of a brane possessing an intrinsic 
curvature term in the action have already been derived and discussed for an AdS-Schwarzschild bulk 
space-time [14]. Cosmological consideration of the DGP model was first discussed in [15] where it was 
shown that in a Minkowski bulk space-time we can obtain self-accelerating solutions. In the original 
DGP model it is known that AD general relativity is not recovered at linearized level. However, some 
authors have shown that at short distances we can recover the AD general relativity in a spherically 
symmetric configuration, see for example [16]. An important observation was made in [17, 18] where 
it was shown that DGP model allows for an embedding of the standard Priedmann cosmology whereby 
the cosmological evolution of the background metric on the brane can entirely be described by the 
standard Friedmann equation plus energy conservation on the brane. This was later extended to 
arbitrary number of transverse dimensions in [19]. For a comprehensive review of the phenomenology 
of DGP cosmology, the reader is referred to [20] . 

Standard cosmology suggests that as the early universe expanded and cooled, it underwent a series 
of symmetry-breaking phase transitions, causing topological defects to form. It is the study of such 
phase transitions that would pave the way for a better understanding of the evolution of the early 
universe, characterized by the existence of a quark-gluon plasma undergoing a phase transition. In 
what follows we focus attention on possible scenarios which might have occurred to allow the phase 
transition mentioned above come to the fore. We generally follow the discussion presented in [21, 22] 
which puts the quark-gluon phase transition in a cosmologically transparent perspective. 

The existence of a phase transition from the quark-gluon plasma phase to hadron gas phase is a 
definite prediction of QCD. However, the phase transition in QCD can be characterized by a truly 
singular behavior of the partition function leading to a first or second order phase transition, and also 
it can be only a crossover with rapid changes in some observables, strongly depends on the values of 
the quark masses. The possibility that the bags of quark-gluon undergoing a phase transition was 
demonstrated for the first time in [23]. Most studies have shown that one may obtain a first, second 
or higher order transition, with the possibility of even no phase transition being pointed out in [24]. 
Recently, lattice QCD calculations for two quark fiavors suggest that QCD makes a smooth crossover 
transition at a temperature of Tc ~ 150 MeV [25]. Such a phase transition could be responsible for the 
formation of relic quark-gluon objects in the early universe which may have survived. In this paper, 
our study on phase transition are based on the ideas proposed in the first reference in [23] where it 
was shown that under certain conditions a gas of extended hadrons could produce phase transitions 
of the first or second order, and also a smooth crossover transition that might be qualitatively similar 
to that of lattice QCD. 
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The cooling down of the color deconfined quark-gluon plasma below the critical temperature, 
believed to be around Tc ~ 150 MeV, makes it energetically favorable to form color confined hadrons 
(mainly pions and a tiny amount of neutrons and protons, since the net baryon number should be 
conserved). However, the new phase does not show up immediately. A first order phase transition 
generally needs some supercooling to overcome the energy expended in forming the surface of the 
bubble and the new hadron phase. When a hadron bubble is nucleated, latent heat is released and a 
spherical shock wave expands into the surrounding supercooled quark-gluon plasma. This causes the 
plasma to reheat, approaching the critical temperature and preventing further nucleation in a region 
passed by one or more shock fronts. Bubble growth is generally described by deflagrations where a 
shock front precedes the actual transition front. The nucleation stops when the whole Universe has 
reheated to Tc- The phase transition corresponding to this phase ends promptly, in about 0.05 /xs, 
during which the cosmic expansion is completely negligible. Afterwards, the hadron bubbles grow at 
the expense of the quark phase and eventually percolate or coalesce. Ultimately, when all quark-gluon 
plasma has been converted into hadrons, neglecting possible quark nugget production, the transition 
ends. The physics of the quark-hadron phase transition and its cosmological implications have been 
extensively discussed in the framework of general relativistic cosmology in [26]- [39]. 

As is well known, the Fricdmann equations in brane- world scenarios differs from that of the 
standard 4£) cosmology in that it results in an increased expansion rate at early times. We expect this 
deviation from the standard AD cosmology to have noticeable effects on the cosmological evolution, 
especially on cosmological phase transitions. In the context of brane-world models, the first order 
phase transition was studied in [40] where it was shown that due to the effects coming from higher 
dimensions, a phase transition requires a higher nucleation rate to complete and baryogenesis and 
particle abundances could be suppressed. Recently, the quark-hadron phase transition was studied in 
a Randall- Sundrum brane-world scenario [21]. Within the framework of phase transition, the authors 
studied the evolution of the relevant cosmological parameters (energy density, temperature, scale 
factor, etc.) of the quark-gluon and hadron phases and the phase transition itself. It would therefore 
be interesting to study such phase transitions in the context of a DGP brane-world scenario with a 
cosmological constant A, hereafter referred to as ADGP, and this is what we intend to do in what 
follows. 

2 Field equations in DGP brane scenario 

We start by writing the action for ADGP brane-worlds as 

S = ^f d''x^n+ [ d^X^(^R-x)+Sm[q^..,^m], (1) 
^ J Bulk J brane \ ^ / 

where the first term is the Einstein-Hilbert action in 5D with the 5-dimensional bulk metric g^^g and 
Ricci scalar TZ. Similarly, the second term is the Einstein-Hilbert action corresponding to the induced 
metric q^i, on the brane, where R is the relevant scalar curvature and A is the cosmological constant 
on the brane, with M4 = (87rG'4)~2 and M5 = (87rG5)^3 being the reduced Planck masses in four 
and five dimensions respectively and 5,71 is the matter action on the brane with matter field ^|Jm■ The 
induced metric q^i, is defined as usual from the bulk metric g^^g by 

Qf.. = Sis!^9^s- (2) 
Variational of the action with respect to g^g yields the field equations 

Gab + {2£^^pS^J-'bG^, - X6^J'^,q^,)5{y) = Tab, (3) 



3 



where T4B is the energy-momentum tensor of the matter content of the 5-dimensional space-time and 
•^DGP ~ 2M^ crossover length scale between the iD and 5D regimes in the DGP brane model. 

Note that the brane is the the hypersurface defined at y = 0. To continue, we follow the approach 
presented in [15] in which it was first shown how one recovers a self-accelerating solution from the 
DGP field equations (3). 

Being interested in cosmological solutions, we take a metric of the form 

= -n\T, y)dT^ + a\T, y)6ijdx'dx^ + b^T, y)dy^ , (4) 

where 7.^ is a maximally symmetric 3-dimensional metric (k = —1,0, 1 will parameterize the spatial 
curvature) and we focus our attention on a spatially flat {k = 0) space-time. Since matter resides on 
the brane and the bulk is taken to be empty, the field equations (3) in the bulk can be written as 
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where a prime denotes differentiation with respect to y and a dot denotes differentiation with respect 

to r. By using a technique first developed in [5, 41], the bulk equations, remarkably, may be solved 
exactly given that the bulk is devoid of matter. Assuming Z2 symmetry across the brane, we find 
that the metric components in equation (4) are given by [15] 



niT,y) = 1 ^ \y\-^ , 

a{r,y) = Oo T \y\do , 
Kr,y) = 1, 



(9) 



where aQ(T) is yet to be determined. Note that ao(T) = a(r, y = 0) represents the usual scale factor of 
the four-dimensional universe. According to equations (9), there are two distinct possible cosmologies 
associated with choice of the sign. 

To obtain cosmological equations we must take the total energy-momentum tensor which includes 
matter and the cosmological constant on the brane. We thus write 

Tfi Ibrane = S{y) dilig{-p,p,p,p,0) . (10) 

In order to determine the scale factor aQ(T), we need to employ the proper boundary condition on 
the brane. This can be done by taking equations (3) and integrating across the brane. Then, the 
boundary conditions require 

87rG4 
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Comparison with the bulk solutions (9) imposes a constraint on the evolution of aQ(T). Such an 
evolution is tantamount to a new set of Friedmann equations [15] 



H'--^H=^-^p{r) + X. (13) 



•^DGP 

and 



p + 3{p + p)H = 0, (14) 

where e = it 1 and we have used the usual Hubble parameter H = a^/aQ. The second of these equations 
is just the usual expression of energy-momentum conservation. The first equation, however, is indeed a 
new Friedmann equation which is a modification of the ordinary four-dimcnsional Friedmann equation 
with a cosmological constant on the brane. As for a value for the brane cosmological constant A, 
several bounds have been suggested. Thus from the big bang nucleosynthesis constraints it follows 
A > 1 MeV^ [42]. A much stronger bound for A arises from null results of submillimeter tests of 
Newtons law, leading to A > 10^ GeV^ [43]. An astrophysical lower limit on A, which is independent 
of the Newton's law and cosmological limits has been derived in [42] with the result A > 5 x 10^ 
MeV^. 

Equations (13) and (14) together with the brane equation of state are sufficient to derive the 
cosmological evolution of the brane metric. Now, assuming p > 0, equation (13) can be rewritten as 

Let us examine equation (13) more closely. The new contribution from the DGP brane- world scenario 
is the appearance of the term ^-^ — H on the left-hand side of the Friedmann equation. The choice for 
the sign of e represents two distinct cosmological phases. Just as gravity is modified at large scales and 
behaves 4£)-like at short distances, so too is the Hubble scale, H(t), which evolves by the conventional 
Friedmann equation at high Hubble scales but is altered substantially as H(t) approaches ^"^p. There 
is therefore the possibility of two distinct cosmological phases [15]. First, the phase in equation (13) 
corresponding to the choice e = — 1 which had already been studied in [14, 44, 45], with transitions 
from ~ /9 to ^ p^. We refer to this phase as the Friedmann- Robertson- Walker (FRW) phase. 
The other cosmological phase corresponds to e = 1. Here, cosmology at early times again behaves 
according to the conventional four-dimensional Friedmann equation, but at late times approaches a 
brane self-accelerating phase. 



3 Quark— hadron phase transition 

The quark-hadron phase transition is a notion fundamental to the study of particle physics, particu- 
larly in the context of lattice gauge theories. However, it has become increasingly relevant and even 
an integral part of any study dealing with the underlying mechanisms responsible for the evolving 
universe at its early stages of formation in which a soup of quarks and gluons interact and undergo 
a first order phase transition to form hadrons. It is therefore essential to have a review of the basic 
ideas before attempting to use the results obtained from such a phase transition and apply them to 
the study of the evolution of the early universe within the context of a DGP brane-world scenario. In 
this regard, a well written and concise review can be found in [21] and the interested reader should 
consult it. Here, it would suffice to mention the results relevant to our study and leave the details of 
the discussion to the said reference. 
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We start from the equation of state of matter in the quark phase which can generally be given in 
the form 

Pq = ^aqT^ + V{T), pg = agT^-V{T), (16) 

where Uq = {Tr^/90)gq, with Qq = 16 + {21/2)Nf + 14.25 = 51.25 and Np = 2 with V{T) being the 
self-interaction potential. For V{T) we adopt the expression [37] 

V{T) = B + ^tT^ - arT"^, (17) 

where B is the bag pressure constant, ax = 77r^/20 and 'yr = with rus the mass of the strange 

quark in the range nig G (60 — 200) MeV. In the case where the temperature effects are ignorablc. 
the equation of state in the quark phase takes the form of the MIT bag model equation of state, 
Pq = {Pq — 4:B)/3. Results obtained in low energy hadron spectroscopy, heavy ion collisions and 
phenomenological fits of light hadron properties give 

5^4 between 100 and 200 MeV [46]. 
Once the hadron phase is reached, one takes the cosmological fluid with energy density ph and 
pressure ph as an ideal gas of massless pions and nucleons obeying the Maxwell-Boltzmann statistics. 
The equation of state can be approximated by 



Ph{T) = ^Ph{T) = a^T\ 

where a^r = (7r^/90)g'/j and gh = 17.25. 

The critical temperature Tc is defined by the condition Pq{Tc) 



(18) 



2{aq + aT - a^) 



Ph{Tc) [47], and is given by 

1/2 



(19) 



If we take iris = 200 MeV and B^/^ = 200 MeV, the transition temperature is of the order Tc ~ 125 
MeV. It is worth mentioning that since the phase transition is assumed to be of first order, all the 
physical quantities exhibit discontinuities across the critical curve. 



4 Dynamical consequences of DGP brane universe during quark— 
hadron phase transition 

We are now in a position to study the phase transition described above. The framework we are working 
in is defined by the DGP brane-world scenario for which the basic equations were derived in section 
2 and therefore the essential physical quantities that should be monitored through the quark-hadron 
phase transition are the energy density p, temperature T and scale factor . These parameters are 
determined by the Friedmann equation (13), conservation equation (14) and the equations of state, 
namely (16), (17) and (18). One should now consider the evolution of a DGP brane-world before, 
during and after the phase transition era. 

Let us consider the era preceding the phase transition for which T > Tc and the universe is in the 
quark phase. Use of equations of state of the quark matter and the conservation of matter on the 
brane, equation (14), leads to 

^ _ a'p _ Sag -a^ f 1 jj. f 

Oq 3aq T 6aqT^' ^ ^ 

Integrating the above equation immediately gives 

«.(T)=cT^exp(-2L_), (21) 
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where c is a constant of integration. 

We may now proceed to obtain an expression describing the evolution of temperature of the DGP 
brane universe in the quark phase by combining equations (20) and (13), giving 



(22) 

where we have denoted 

Ao = l-^, (23) 

Equation (22) may now be solved numerically and the result is presented in Figure 1 which shows the 
behavior of temperature as a function of the cosmic time r in a DGP brane-world filled with quark 
matter for different values of A. 

4.1 Constant self-interaction potential (y{T) = B) 

One may gain an analytical understanding of the evolution of cosmological quark matter in our DGP 
brane-world by taking the negative branch (e = —1) in equation (13). Let us concentrate on the 
phase transition era, namely <C ^^Qp and take the simple case in which temperature corrections 
can be neglected in the self-interaction potential V . Then V = B = cons, and equation of state of 
the quark matter is given by that of the bag model, pq = {pq — AB) /3. Equation (14) may then be 
integrated to give the scale factor on the brane as a function of temperature 

aoiT) = ^, (24) 

where c is a constant of integration. 

Since in the early universe the density of matter is very high and the Hubble radius {H~^) is 
small with respect to the crossover length scale (^^gp) ^ ^dgp ' equation (13) for e = — 1 

leads to i.e. ~ p^ which is the specific character of brane world models. The evolution of the 
quark phase of the DGP brane-world is now described by the equation 

<=^p + X. (25) 

Og 3 

Using equations (14) and (25), the time dependence of temperature can be obtained from equation 

dT -87rG4 



dr 3 

This equation has the following general solution 

(_8^5 + A)V4 

[SirG^aq - Cexp(^^^f^ - 4A)r 



{3aqT^ + BT) - XT. (26) 



= ^ ' — T74- (27) 



The integration constant C can be obtained by relating it to the temperature T'|,-=t-o = of the 
quark matter at the time tq 



SirG^aqT^ + ^B + A 
exp(5^^ - 4A)to 
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Figure 1: The behavior of T{t) as a function of time (r) for the positive branch (e = 1) and different values of A: 
A = 8 X 10* MeV* (solid curve), A = 8 x 10'' MeV** (dashed curve), A = 2 x 10^° MeV* (dotted curve) and A = 1 x 10" 
MeV* (dotted-dashed curve). We have taken B^^'^ = 200 MeV , G4 ~ lO"^'' MeV"^ and ^^gp ~ lO^^cm [12]. 

4.2 Formation of hadrons 

During the phase transition, the temperature and pressure are constant and quantities hke the entropy 
S = sa^ and enthalpy W = (p + p) are conserved. Following [47, 21], we replace p (r) by /i(t), the 
volume fraction of matter in the hadron phase, by defining 

p (r) = p^h{T) +Pq[1- Ht)] =Pq[1 + nh{T)] , (29) 

where n = (/o^ — /Jq) /pg. The beginning of the phase transition is characterized by /i(tc) = where 
Tc is the time representing the beginning of the phase transition and p (tc) = Pg , while the end of the 
transition is characterized by h (rh) = 1 with Th being the time signaling the end and corresponding 
to p{Th) = p^. For T > Tfi the universe enters into the hadronic phase. Equation (14) now gives 

< _ 1 {Ph-Pq)^ _ 1 rh ^^^^ 



«o PQ+Pc+ [Ph - Pq)^ 3 1 + r/i' 

where we have denoted r = (p^ ~ Pq) I {Pq +Pc)- From the above equation we can obtain the 
relation between the scale factor on the brane and the hadron fraction /i(r) 

a,{T)=a,{T,)[l + rh{T)]-'/\ (31) 

where the initial condition h (tc) = has been used. Now, using equations (13) and (31) we obtain 
the time evolution of the matter fraction in the hadronic phase 




+ X + ^Pg[l + nh{r)]j. (32) 

Figure 2 shows variation of the hadron fraction /i(t) as a function of r for different values of A 

It is well known in brane-world models that a quadratic term in the energy density appears in 
the Friedmann equation and that in the DGP scenario such a term can be obtained by taking the 
negative branch (e = — 1) when H <^ ^dgp' describing the effects of the extra dimensions. Thus, at 
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h(r) 




Figure 2: The behavior of /i(r) as a function of time (r) for the positive branch (e — 1) and different values of A: 
A = 8 X 10* MeV (solid curve), A = 8 x 10^ MeV* (dashed curve), A = 2 x 10^° MeV (dotted curve) and A = 1 x lO" 
MeV (dotted-dashed curve). We have taken B^/^ = 200 MeV , G4 - lO"** MeV'^ and ^dgp ~ lO^^cm. 



the beginning the quadratic term dominates in the evolution of the universe so that the evolution of 
the hadronic fraction during the phase transition is given by the equation 

dh _ (1 + rh) 
dr r 
with the general solution given by 



-[87rG4/9Q(l + n/i) + A], 



exp[(87rG4PQ + A)t + «^C'] - {SttG^Pq + A) exp 



+ rC'{8.G,p, + A)' 



nSirGiPg exp 



T + rC'{8TTGiPQ + A) - r exp[(8^G4PQ + A)r + 



The integration constant C can be obtained by relating it to h (tc) = at the time tq 



SttG^PqUTq + rln 



^(87rG4Pg+A)To 
87rG4p„+A 



r(87rG4P„n — 87rG4p„r — Ar) 



(33) 



(34) 



(35) 



From equation (31) we see that the evolution of the scale factor on the brane associated with the 
hadronic fraction during the phase transition can be obtained. We also note that the end of the phase 
transition era corresponds to /i(r) = 1, thus one obtains the time at which the phase transition 
ends. 



4.3 Pure hadronic era 

At the end of the phase transition the scale factor of the universe has the value 

a,{Th)=a,{T,){r + l)^ . (36) 

Also, the energy density of the pure hadronic matter after the phase transition is ph = 3ph = Sa-^T"^. 
The conservation equation on the brane (14) gives a^^T) = Ug {th) Tc/T. The temperature dependence 
of the DGP brane universe in the hadronic phase is governed by the equation 
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Figure 3: The behavior of T{t) as a function of time (r) for the positive branch (e = 1) and different values of A: 
A = 8 X 10** MeV* (solid curve), A = 8 x 10^ MeV** (dashed curve), A = 2 x 10^° MeV* (dotted curve) and A = 1 x lO" 
MeV* (dotted-dashed curve). We have taken B^/* = 200 MeV , G4 ~ 10"** MeV'^ and ^^gp ~ lO^^cm. 

Variation of temperature of the hadronic fluid filled ADGP brane universe as a function of r for 
different values of A is represented in figure 3. 

5 Bubble nucleation 

In this section, we consider the process of formation and evolution of microscopic quark nuclei in the 
cosmological fluid on the ADGP brane in the context of bubble nucleation theory. The main target 
of nucleation theory is to calculate the probability that a bubble or droplet of the A phase appears 
in a system initially in the B phase near the critical temperature [48]. 

In the previous section we discussed the mechanisms of the phase transition and thus for the sake 
of completeness, we should consider the phase transition in the context of bubble nucleation theory 
in this section. Phase transition can be described by an effective nucleation theory in the small 
supercooling regime. As the temperature decreases, there is a probability that a droplet of hadrons 
is nucleated from the quark plasma. The nucleation probability density is given by [47] 



where Pg, are dimensionless constants and T = T/Tc- When a droplet of hadrons is formed, it 
starts to expand explosively [49], with a velocity v^^ smaller than the sound speed. At the same time, 
a much quicker shock wave is generated. During the period over which the temperature decreases, 
a large number of bubbles are created until the shock waves collide and re-heat the plasma to the 
critical temperature. The fraction of the volume which at a time r is turned into the hadronic phase 
in the small supercooling scenario is given by [21] 



where Tj is the time at which the critical temperature is reached, and the nucleation process started. 

Let us consider the effects of the DGP brane scenario in bubble nucleation in the early universe. 
At early times, the quadratic term of the energy density in Friedmann equation was dominant. Thus, 




(38) 




(39) 
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Figure 4: The behavior of f{T) as a function of the normaHzed temperature (T) for the different values of A: A = 8 x 10* 
MeV* (sohd curve), A = 8 x 10^ MeV* (dashed curve), A = 2 x 10^° MeV^ (dotted curve) and A = 1 x lO" MeV* 
(dotted-dashed curve). We have taken B^/** = 200 MeV , G4 ~ lO-"*"* 



MeV"^ and l^^^ ~ lO^^cm. 



using equations (25) and (26), we can express this integral in terms of the normahzed temperature 
and calculate the it by a saddle point approximation. The integral (39) may now be approximated as 



/(T) ^ X{f) 



exp 



\1~{T+X{T))* 



(8^)2(5 + A)3 (T + X(f )) [3a^{T + X(r))4 + b] 



X < 



log 



3a,(r + X(f))4 + 6(^±£(Sy 
3a^if + X{f))^ + b 



where b = and X{T) is given by 



(40) 



X{f) 



'l-f^f{6agf^ + b) 



(41) 



In figure 4 we have plotted f{T) against T for some values of A with Vsh = 10"'^ and wq = po = 1- 
The fraction of hadronic matter stays very close to zero until the supercooling temperature is between 
T = 0.9856 and T = 0.9853, thus near T^., it suddenly jumps to 1. Here we may deduce the result 
that at a certain temperature below the critical value, an enormous amount of hadronic bubbles 
are nucleated almost everywhere in the plasma, which grow explosively to transform all the plasma 
into hadrons. This result is similar to what happens in standard cosmology [47] where the small 
supercooling phase transition is also very rapid with respect to the simple first order phase transition, 
at temperature T ~ 0.9856rc, which is very similar to the one we have obtained in the context of 
DGP brane world scenario. Another remarkable feature that is present in the supercooling scenario 
is the dependence of the transition temperature on the cosmological constant, A, in the ADGP model. 



6 Conclusions 

In this paper, we have discussed the quark-hadron phase transition in a DGP brane cosmological 
scenario in which our universe is a three-brane embedded in a 5-dimensional bulk space-time with an 
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intrinsic Ricci curvature scalar and a cosmological constant being included in the brane action, within 
an effective model of QCD. We have studied the evolution of the physical quantities, relevant to the 
physical description of the early universe; the energy density, temperature and scale factor, before, 
during, and after the phase transition. In particular, due to the high energy density in the early 
universe, we have studied the specific case of neglecting the terms linearly proportional to the energy 
density with respect to the quadratic terms in the Priedmann equation when the negative branch 
is considered. By assuming that the phase transition may be described by an effective nucleation 
theory, we also studied the case where the universe evolved through a mixed phase with a small initial 
supercooling and monotonically growing hadronic bubbles. 
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